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We calculate high-harmonic generation (HHG) by intense infrared lasers in atoms and molecules 
with the inclusion of macroscopic propagation of the harmonics in the gas medium. We show that 
the observed experimental spectra can be accurately reproduced theoretically despite that HHG 
spectra are sensitive to the experimental conditions. We further demonstrate that the simulated (or 
experimental) HHG spectra can be factored out as a product of a "macroscopic wave packet" and the 
photo-recombination transition dipole moment where the former depends on the laser properties and 
the experimental conditions, while the latter is the property of the target only. The factorization 
makes it possible to extract target structure from experimental HHG spectra, and for ultrafast 
dynamic imaging of transient molecules. 

PACS numbers: 42.65.Ky,31.70.Hq,33.80.Eh 



High-harmonic generation (HHG) has been employed 
to probe electronic structure of molecules on an ultra- 
fast time scale in recent years [iHSf- When molecules are 
placed in an intense laser field, electrons that are removed 
earlier may be driven back to recoUide with the parent 
ion. HHG occurs when the returning electrons recombine 
with the parent ion with the emission of high-energy pho- 
tons as in an inverse photoionization (PI) process. Since 
PI is a sensitive tool for probing electronic structure of 
molecules, HHG may serve likewise, but with the advan- 
tage of ultrafast temporal resolution as well as covering a 
coherent broad spectral range from XUV to soft-X-rays. 
Experimentally, however, HHG is generated from all the 
molecules in the interaction region. The radiations from 
them co-propagate with the fundamental infrared (IR) 
beam nonlinearly. To extract structure information of 
individual molecules, e.g., the amplitude and phase of 
PI transition dipole from the measured HHG, the prop- 
agation effect in the medium should be investigated. For 
molecular targets, this has not been done so far. Instead, 
it was often assumed that HHG was measured under the 
perfect phase-matching conditions and that the observed 
harmonics were directly proportional to the harmonics 
from a single molecule. While such assumptions may 
be adequate for explaining many experimental observa- 
tions qualitatively, such as the dependence of HHG on 
molecular alignment and on symmetry of the molecular 
orbital, the two-center interference (j, i5j , and multiple- 
orbital contributions to HHG they are inadequate 
if accurate structure information of individual molecules 
are to be extracted from the observed HHG spectra. 

The effect of macroscopic propagation on the observed 



HHG spectra for atoms has been investigated extensively 
in the past two decades in connection with the genera- 
tion of attosecond pulses. However, we are not aware of 
any direct comparison between experimental HHG spec- 
tra and theoretical simulations over an extended spec- 
tral region. The Maxwell's equations that govern the 
propagation of the fundamental driving IR field and the 
generated harmonics are well established. To carry out 
such propagation calculations, accurate induced atomic 
dipoles generated by lasers for hundreds of peak intensi- 
ties should be calculated which serve as the source term 
of the harmonics. These induced dipoles are often cal- 
culated using the strong-field approximation (SFA), or 
the so-called Lewenstein model [7|]. The SFA does not 
describe the laser-atom interactions accurately, thus the 
results from the propagation can only be used to qualita- 
tively interpret the experiments. While accurate induced 
dipoles can be obtained from solving the time-dependent 
Schrodinger equation (TDSE), the calculation is rather 
time consumingand was rarely attempted except for very 
few occasions [8|. Thus after two decades, our under- 
standing of experimental HHG data is still mostly at the 
qualitative level. 

In this Letter, we show that such limitations have been 
removed. We generate HHG spectra theoretically under 
experimental conditions and the results are compared di- 
rectly to the observed data. The simulated spectra agree 
well with the measured one, over a broad photon energy 
region. The experiments were taken in a geometry as de- 
picted in Fig.[TJa), using IR laser pulses with wavelengths 
of 1200 nm and 1360 nm, respectively, generated in a 
high-energy optical parametric amplifier (HE-TOPAS). 
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The TOPAS was pumped by the output of a multipass 
femtosecond amplifier system (8 mJ, 32 fs, 50Hz). High- 
harmonic radiation generated from the gas jet was al- 
lowed to propagate and a slit was placed downstream. 
After the slit, a concave grating dispersed the harmonics 
which were then detected with a CCD detector. The gas 
jet was formed from a supersonic expansion of Ar or N2 
at a stagnation pressure of 3 bars. Experiments were per- 
formed on either isotropic or aligned N2 molecules. The 
experimental spectra were corrected for the response of 
the grating and detector. 

Fig. [TJb) shows the HHG spectra of Ar generated by 
a 1200 nm laser. The horizontal axis is the photon en- 
ergy and the vertical axis is the transverse spatial dimen- 
sion. The upper frame is from the measurement, while 
the bottom frame is from the simulation. The two spec- 
tra are normalized to each other at harmonic 75, or at 
photon energy of 77 eV. There is a general agreement be- 
tween the two spectra. The "up-down" asymmetry in the 
experimental HHG spectra is due to asymmetry in the 
laser beam profile. The faint features near 50 eV are the 
"famous" Cooper minimum in Ar 9], observed in pho - 
toionization, as well as in earlier HHG spectra 13, HI- 
The harmonic yields integrated over the vertical dimen- 
sion are compared in the upper half of Fig. [UJc). The 
lower half shows the HHG spectra taken with the 1360 
nm laser. In both cases, we can see very good agreement 
between theory and experiment over the 30-90 eV region 
covered. Experimentally, the gas jet is 0.5 mm long and 
placed 3 mm after the focus. A vertical slit with a di- 
ameter of 100 /im is placed 24 cm after the gas jet. For 
the 1200 (1360) nm the beam waist at the laser focus 
is 47.5 (52.5) ^m, and the pulse duration is ~40 (^^50) 
fs. To achieve best overall agreement, in the simulation 
the peak intensity and gas pressure for each wavelength 
are adjusted till best overall fit in the data are achieved. 
Thus for 1200 nm laser, the peak intensity for experiment 
(theory) is 1.6 (1.5) x lO^-^W/cm^, gas pressure is 28 (84) 
Torr. For the 1360 nm laser, the corresponding param- 
eters are 1.25 (1.15) x lO^-^W/cm^, and 28 (56) Torr, re- 
spectively. By using a higher pressure in the simulation, 
we find that the higher harmonics become sharper, as 
in experiments. The (normalized) envelope of the har- 
monic spectra, however, does not depend much on the 
gas pressure, see Fig. [S^b) below. 

In the theoretical simulation, we first obtain single- 
atom induced dipole using the quantitative rescattering 
(QRS) theory p^ - [l3 |. The Ar is treated in the single- 
active electron approximation using model potential pro- 
posed by Muller [l5j. The resulting induced dipoles for 
different peak intensities are then fed into the Maxwell's 
equations. The propagation equations for the fundamen- 
tal field and the harmonics are the standard ones [16t - 
[l8l |. For Ar target, we include dispersion, absorption, 
Kerr and plasma effects on the fundamental field in the 
medium. For the harmonics, only the dispersion and ab- 
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FIG. 1: (Color online) (a) Typical configuration for measur- 
ing HHG in the far field, (b) HHG spectra of Ar generated by 
1200 nm laser. Upper frame: experiment; lower frame: the- 
ory, (c) Comparison of theoretical (green curves) and exper- 
imental (red curves) HHG yields integrated over the vertical 
dimension for 1200 nm (upper curves) and 1360 nm (lower 
curves) lasers. Laser parameters are given in the text. 



sorption are included. The harmonic yields emitted at 
the exit face of the gas jet (near field) are propagated to 
the far field where the harmonics are measured. They are 
obtained from the near-field harmonics through a Han- 



kel transform [19|, |20| . We assume the laser beam in the 
entrance of gas jet has the Gaussian shape. 

A careful examination of Figs.[IJb) and (c) reveals that 
there are still small discrepancies between the experimen- 
tal data and the simulation. The harmonic width (or 
harmonic chirp) is narrower from the theory than from 
the experiment. Harmonic chirp is a direct consequence 
of temporal variation of laser intensity. The harmonic 
width is mainly influenced by the pulse duration, pres- 
sure, and laser intensity [2l| - |23| . The width decreases 
with increasing pulse duration, and with decreasing gas 
pressure. Other experimental factors like use of the slit 
and location of the detector also can affect the HHG spec- 
tra. 

High-order harmonics from molecules by 1200 nm 
lasers have been reported for aligned and randomly dis- 
tributed N2 and CO2 recently [2^. Here we report our 
simulated results for N2, at the two peak laser intensities, 
0.9 and l.lxlO^'*W/cm^, reported in [2J]. To achieve 
good agreement in the cutoff positions, the two inten- 
sities used in the theory are 0.78 and 0.9xlO^^W/cm^ 
instead, respectively. Since the experiment was carried 
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out at low laser intensity and low gas pressure, the har- 
monics are propagated without absorption and disper- 
sion effects from the medium, and the fundamental laser 
field is not modified through the medium [2^. In the 
theoretical simulation, we first obtain induced di pole s of 
fixed-in-space molecules using QRS theory jl^, 14 1 for 
different laser peak intensities. The induced dipoles are 
averaged coherently according to the alignment distribu- 
tion and then fed into the Maxwell's equations. Fig. [2] 
shows the good overall agreement between the measured 
and the simulated spectra, for both randomly distributed 
and aligned N2. By examining the experimental HHG 
spectra more carefully, they reveal a shallow minimum at 
38±2eV (low intensity) and at 41±2eV (high intensity) 
for both aligned and unaligned molecules. The theory 
also predicts a minimum: for unaligned molecules, the 
minimum is at ~39eV for low intensity and ^40eV for 
high intensity. For aligned molecules, the minimum is at 
~42eV for low intensity and ~44eV for high intensity. In 
the experiment, the degree of alignment was estimated 
to be (cos^ 0}=O.6-O.65. In the simulation, an alignment 
distribution of cos** 9 is used. Note that only HOMO is in- 
cluded in the calculation. We believe that this is the first 
time that HHG spectra from molecules have been cal- 
culated including the propagation effect in the medium 
and the simulated results have been compared directly 
to the measured spectra. In the future, HHG spectra 
taken at different alignment angles should be compared 
together. Such comparison would help to identify factors 
that contribute to the remaining discrepancies between 
experiment and simulation. 

The macroscopic HHG spectra can be expressed as [2I] 



(1) 



where W{uj) (the complex amplitude) is called "Macro- 
scopic wave packet" (MWP), d(aj) is PI transition dipole 
moment for the atom. For the molecule d{ijj) is taken 
to be the coherently averaged PI transition dipole mo- 
ment df'''3{uj) = /J" N{e)ip{9)d{9,uj) siii9d9, where N{e) 
is the alignment-dependent ionization probability, p{9) 
is the alignment distribution, and d{9, w) is the parallel 
component of the alignment-dependent transition dipole 
moment (26j . The polarization of pump laser is assumed 
parallel to probe laser. For unaligned molecules, p{9) 
is a constant. Actually, MWP has the clear physical 
meaning. It can be considered as the collective effect 
of microscopic wave packets for the returning electrons 
I3I 14 1, which is governed by Maxwell's equations. In 
other words, the laser and macroscopic medium effects 
are all combined into MWP. 

The validity of Eq. ([1]) has been checked in Jin et al. 
[isj when both the laser intensity and the gas pressure 
are low. The correctness of this relation has been as- 
sumed in Itatani et al. [l[ by comparing Ar with N2, and 
in Levesque et al. for rare gas atoms. However, this 
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FIG. 2: (Color online) Comparison of HHG spectra from the- 
ory (green curves) and experimental results (red curves) of 
Ref. [ill, (a) for randomly distributed N2 and (b) for N2 
aligned along laser polarization direction. The laser intensi- 
ties are indicated where /o=10^^ W/cm^. See text for addi- 
tional laser parameters. 



relation has not been carefully checked for different fo- 
cusing conditions and laser parameters. Theoretically, we 
have checked the validity of Eq. ([TJ carefully. We carried 
out macroscopic propagation calculation of HHG using 
single-atom (single-molecule) induced dipole obtained by 
QRS and by SEA. We have been able to show that the 
MWP obtained from the two calculations agree rather 
well, irrespective of laser parameters or the focusing con- 
ditions. In other words, the medium propagation only 
affects HHG through its modifications on the MWP. In 
this way, to study propagation effect on HHG, we can 
just study how the MWP (only the amplitude is consid- 
ered below) depends on the lasers and the experimental 
conditions. 

In Fig. [3l^a) we show the dependence of MWP on the 
position of the Ar gas jet with respect to the laser fo- 
cus. The laser intensity is 1.6x 10^'*W/cm^, and the gas 
pressure is 56 Torr. For easy visualization we show the 
smooth envelope of |VF(a;)|. The three curves are for the 
gas jet at z=-3 mm (laser focus before gas jet), z=0 (at) , 
-1-3 (after). It is generally known that HHG achieves best 
phase-matching if the gas jet is placed behind the laser 
focus where the dipole phase from the harmonic can be 
partially canceled by the Gouy phase. Thus among the 
three curves, the "after" curve is the flattest one. For the 
"before" focus, the MWP varies most as the photon en- 
ergy is changed, and the phase (not shown) varies widely 
from order to order — reflecting poor phase-matching for 
this geometry. 

In Fig. El^b) we compare the MWP derived from chang- 
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FIG. 3: (Color online) (a) Dependence of macroscopic wave 
packet (MWP) on the position of the Ar gas jet with 

respect to laser focus; (b) Effect of gas pressure on MWP. 
The curves are renormalized such that they should fall on the 
same curve if the perfect phase-matching condition is fulfilled, 
(c) The MWP for the two lasers, and the magnitude of the 
PI transition dipole moment of Ar. (d) Same as in (c) but for 
N2. The MWP is for two different laser intensities, and the 
averaged PI transition dipoles are for isotropic and aligned 
molecules. The laser intensities are indicated where /o = 10^* 
W/cm^. See text. 



ing the Ar gas pressure for the "after" focusing condi- 
tion. The MWP has been normalized by the ratio of 
the pressure. Under perfect phase-matching condition, 
the MWP (the amplitude) should be propor- 

tional to the pressure [2^, [2^. The three curves are 
on top of each other from 45-75 eV, indicating good 
phase-matching in this energy region, but differ some- 
what at lower and higher energies, indicating good phase- 
matching condition is not fulfilled. This demonstrates 
that phase-matching condition cannot be achieved for all 
the harmonics in a given experiment. 

According to Eq. ([1]), the minimum in the HHG spectra 
can occur for different reasons. In Fig. ^c), the MWP 
derived from Ar target using 1200 nm and 1360 nm lasers 
are shown. The two MWP's are quite similar but near 50 
eV, they have slight different slopes. On the other hand, 
the PI transition dipole reveals a clear but broad Cooper 
minimum near 50 eV. Thus the broad minimum in the 
HHG spectra shown in Fig. [TJc) is due to the minimum 
in the PI transition dipole. To pin down the position of 
the "real" minimum, on the other hand, is not as easy 
since the minimum position can be modified somewhat 
by the MWP. 

Similar analysis can be carried out on the HHG spec- 
tra of N2 shown in Fig. [2] The averaged PI transition 
dipole indeed shows a rapid drop near 40 eV, which is 
due to the presence of a shape resonance [s^] of N2 in 
the lower energy. The rapid drop is more pronounced for 
aligned molecules than for random ones, see Fig. [2Id). 



For the MWP, under the same laser intensity, we have 
checked that they are the same for randomly distributed 
and aligned molecules. Thus it explains why the HHG 
from single-molecule response can be used to interpret 
how the intensity of each harmonic changes with pump- 
probe time delay in Le et al. 3l|. However, the MWP 



changes more rapidly with laser intensity, especially for 
the longer wavelength laser used here. We note that the 
two MWP's in Fig. |3Jd) have somewhat different slopes 
near 40 eV. The multiplication of the MWP and the 
PI transition dipole results in a weak minimum in the 
observed HHG spectra. The minimum would be more 
clearly seen if the molecules were better aligned. From Le 
et al. [l3| . the minimum in PI transition dipole changes 
rapidly with the alignment angle and the eff'ect is severely 
averaged out when molecules are not well aligned. We 
further mention that the MWP in Figs, ^c) and (d) are 
rather different. They are due to the large difference in 
the laser peak intensities used [2^ . In the future, it is de- 
sirable that predictions such as those in Fig. [3] be checked 
experimentally. 

In summary, we demonstrated that experimental HHG 
spectra can now be accurately reproduced from ab initio 
calculations. The theory starts with the calculation of 
laser-induced dipole from single atom or molecule using 
the recently developed quantitative rescattering (QRS) 
theory 13|, [lj| . The propagation effect of the fundamen- 
tal field and the harmonics in the medium is incorporated 
by solving the Maxwell's equations. We further showed 
that the simulated (and experimental) HHG spectra can 
be expressed as the product of a "macroscopic wave 
packet" (MWP) and the photo-recombination transition 
dipole moment. The latter is a property of the target, 
and is independent of the lasers, nor of the propagation 
effect. This factorization makes it possible to extract 
target structure information from the experimental HHG 
spectra. It provides the needed theoretical basis for us- 
ing HHG as ultrafast probes of excited molecules, such as 
those demonstrated recently Clearly, this work also 
opens up opportunities for the quantitative studies of the 
phases of HHG which are fundamental to the generation 
of attosecond pulses. 
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